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Two real vector fields are revealed as spin connections of the spinor field, which is introduced as 
a representation of the local Lorentz group by Dirac spinors. One of these fields is identified as the 
Maxwell field. Another one is the axial field, which is sufficient to describe parity nonconservation 
effects in atoms and it is sensitive to the second internal degree of freedom of the Dirac spinor field. 
The polarization associated with this degree of freedom is a distinctive characteristic of the Dirac 
spinor field in spatially closed localized states. The short-distance behavior of the axial field can lead 
to a Dirac Hamiltonian that is not self-adjoint and trigger a "falling onto a centre" phenomenon. 
In the second part of this work, the axial field will be linked to the gravity. 

PACS numbers: 03.30.-|-p, 03.65.-w, 03.70.-|-k, 04.40.-b,12.10.-g,98.80.-k 



I. INTRODUCTION 

In the late 1920's, V. Fock formulated a vigilant ap- 
proach P] to the problem of "minimal interaction" be- 
tween the spinor matter and the Maxwell field, which 
was not based on Weyl's ad hoc declaration of gauge in- 
variance 2] . Fock derived the Maxwell field as a part of 
the spin connection of the Dirac field. The more confi- 
dence we gain that spinor fields should be thought of as 
the fundamental form of matter, the more rationale we 
have to explore an immense potential of a so consistent, 
economic and elegant method. There is no doubt that 
the spinor fields are the most important objects in the 
physics of stable matter (the Pauli principle is decisive 
for its actual existence). Moreover, the two-component 
spinors represent the light cone of special relativity in 
a most straightforward way, for their currents are light- 
like. As a result, the dynamics of spinors dominate those 
physical phenomena where sharp space-time localization 
occurs. 

In this paper I show that, following the Fock method, 
one can derive the existence of a second real vector field, 
which is minimally coupled to the axial vector current 
of the Dirac spinor field. From a theoretical perspective, 
this field seems to be just overlooked, and (in I920's) 
there was no experimental evidence of its existence. For 
the lack of an established term, I shall refer to it as an 
axial field. This field is not gradient invariant and thus 
should be massive. At the first glance, it very much re- 
sembles the neutral vector field of the electro-weak the- 
ory. The axial field resolves the second internal degree of 
freedom of the Dirac spinor field, the existence of which 
was pointed out by Fock |^ . The presence of the axial 
field in the Dirac equation reasonably describes the (well 
known by now) parity non-conservation (PNC) phenom- 
ena in atomic physics. The most distinctive physical fea- 
ture of this axial field is that its Lorentz force pulls the 
left and right components of a 4-spinor in opposite direc- 
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tions. This causes an additional polarization of the Dirac 
field and additional types of spinor modes (with quantum 
numbers that correspond to compact objects) become 
possible. Like the electromagnetic field, the axial field 
emerges as one of the ingredients of parallel transport of 
Dirac spinors; it is a kinematic effect stemming from the 
complex nature of the spinor field and a large number 
of its polarization degrees of freedom. All immediately 
anticipated physical effects of the axial field seem to be 
known or expected. Surprising is the fact that this field 
is derived at the most basic level of Lorentz invariance 
and that it was not motivated by any phenomenological 
input, e.g., a symmetry observed in particular processes. 

A major focus of this first paper is on the deriva- 
tion and classical aspects of the axial field. This in- 
cludes a detailed analysis of parallel transport of the 
Dirac spinors, simultaneous derivation of the MaxwcU- 
and axial- fields as the components of the spin connec- 
tions, and an analysis of the effect of the external axial 
field on the spherically-symmetric states (hydrogen-like 
atom). The fact that the axial field makes the Dirac 
equation extremely singular is left aside in this paper. 
This fact is merely demonstrated, which then leads to 
more subtle aspects of the axial field, including its role 
in auto-localization of the Dirac field and its relation to 
gravity, which are discussed in the second paper 0. 



II. PARALLEL TRANSPORT OF THE DIRAC 
FIELD. 

To define the derivative of a vector or tensor field one 
has to simultaneously compare many components at two 
neighboring points. Only for vectors given in Cartesian 
coordinates does the notion of parallel transport have no 
ambiguity. One has to account for a continuous rotation 
of the local coordinate hedgehogs along a displacement 
path and include a connection into the covariant deriva- 
tive of a vector even if the same fiat space-time is parame- 
terized by oblique coordinates . The covariant derivative 
of a vector can then be derived in a relatively simple way 
because rotation of a vector at a given point follows the 
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rotation of the local coordinate axes. There is no equally 
simple rule for spinors. 

A. Tetrad formalism 

The tetrad formalism ^ directly inherits the princi- 
ple of local equivalence of the theory of relativity. The 
four components of a usual (coordinate) vector are 
replaced by four coordinate scalars Va, which behave as 
the vectors with respect to Lorentz rotations (x) of lo- 
cal coordinate axes, = eJ^Va, Va{x) A^{x)Vb{x) . 
The four tetrad vectors are the Lorentz vectors and 
coordinate vectors at the same time. The curvilinear 
metric g^^{x) is connected with the local metric gat of 
Minkowski space by 

gf^A^) - gabe''f,{x)e\{x), (2.1) 
gab ^ gM-(a;)ea^(a;)e''Ja;) = diag[l, -1, -1, -1] , 

and the latter is the same at all space-time points. The 
coordinate and the tetrad components of the covariant 
derivatives are ^\V^ = d\Vf_i — T'^xV^, and DaVb = 
e^(VAV^)eJ^ = daVij—ujbcaV^, respectively. In these equa- 
tions, rj^^ are the Christoffel symbols, da = ej^c?^ is the 
derivative in direction a, and LUabc are the Ricci rotation 
coefficients, ujbca = (V^e|^)ec,ye^ = -ujcba ■ 

It is useful to keep in mind that the absolute differ- 
ential, DVa, of a vector Va is the principal linear part 
of the vector increment with respect to its change in the 
course of parallel transport along the same infinitesimal 
path. (A bonus for this definition is that DVa is also 
a vector.) Therefore, the parallel transport just means 
that DVa = 0, and the same definition should be adopted 
for spinors as long as they can represent vector observ- 
ables. As a matter of fact, the tetrad representation most 
adequately reflects the local nature of vector fields in rel- 
ativistic field theory. The Lorentz transformation and 
the definition of parallel transport for spinors are far less 
obvious, because their components are not directly con- 
nected with the vectors of the coordinate axes. Therefore, 
spinors should always be treated as coordinate scalars. 

All observables associated with the spinor fields are 
Lorentz tensors, which are bilinear forms built with the 
aid of Dirac matrices. Within the tetrad formalism, these 
matrices should be treated as coordinate scalars: They 
are the same at each local Lorentz frame at each world 
point. We use an old-fashioned convention that corre- 
sponds to the spinor representation. The conjugated 
spinor is t/)"*". The sixteen Dirac matrices can be con- 
veniently arranged as the products like paCb, a,b ~ 
0,1,2,3. The basic matrices pi and cr^, i = 1,2,3, 
were introduced by Dirac 0; we reserve various nota- 
tion, po = ao = 1 , for the same unit matrix. The other 
notation for the Dirac matrices are: a" — (q;°, a*)(with 
a° = 1, at = 7°7' = PsCTj), pi = {3 = 7°, p2 = -il^l^, 
ps = —7^. The 4x4 matrices a and p satisfy the 
same commutation relations as the Pauli matrices, and 



all matrices a commute with all matrices p : (TiCFk = 

5ik + iiiklC^k, PaPb = Sab + i^abcPc , f^iPa — PaC^i = 0. 

Three 4x4 matrices ai have a block form with the 
standard 2x2 Pauli matrices r^, 

along the block-diagonal, i.e. ai = 1 (8) t^. Matrices p are 
the same products in inverse order, pa = (g) 1, 

"•"(or.)- = ( i ) ■ 

«=(.°x-'„^)--(;-°0- 



B. Parallel transport of Dirac spinors and spin 
connections 

One can look at the spinors as the fields and differen- 
tiate them only after having stated the rules for parallel 
transport. To derive these rules we must rely on various 
tensors, which can be arranged as a matrix, ip'^ paa^ip- It 
is most natural to require the vector current, 

ja = l/'^aa'0 = 'i/'^(l,P30-i)V' = ^lai> , (2.4) 

to be a Lorentz vector. Indeed, its temporal component is 
the unit quantum mechanical operator, which thus com- 
mutes with all other operators, including the Hamilto- 
nian, and by all means, represents a unique conserved 
quantity. This operator is allied with the time-like tetrad 
vector that designates the direction of Hamiltonian evolu- 
tion and thus defines all observables and, ultimately, the 
physical content of a theory. The vector ja{x) is trans- 
formed as ja{x) — + Aa{x)jb{x) uudcr the local Lorentz 
rotation, and its variation under the displacement dx^ is 
Sjf^ — V^^ji^dx^. The tetrad components of this vector 
change by 

Sja = t^acbfds'' = UJacbi'^Oi'll'ds^ , (2.5) 

when this vector is transported at ds°. Let matrix Fq 
(the spin connection) define the change of the spinor com- 
ponents in the course of the same infinitesimal displace- 
ment, S^l) = Tatpds"', Sip'^ = ^/'"'"r+ds". This gives yet 
another expression for Sja, 

Sja = V+ (T+aa + aarb)^ds'' . (2.6) 

This must be the same as Eq. (|2.5ll . Hence, the equation 
that defines Fa is 

F+tta + UaTb = UJacbO." ■ (2.7) 

The matrices {1, Pa, Ci, paCi} form a complete set, and 
all of them, except for the unit matrix, are traceless; 
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the square of each of these matrices is a unit matrix. 
Therefore, the spin connection can be expanded as 

Ti=ai+J2 r[pa + J2 + E fLp-'^r (2-8) 



This results in a system of equations for the coefficients 
that has the following solution, 



ab 



-ieAb, r3 b ^ -ig^b, 



(2.9) 



where Ab and H;, are the tetrad components of two real 
vector fields. The origin of these fields is very sim- 
ple: Eqs. (|2.7() yield only two relations, Re[a^] = and 



Mrl] 



0, for the complex coefficients and rg, which 



leave their imaginary parts arbitrary functions. Thus, 



1 i 

Vb{x) = -ieAb{x) - igps'i^bix) - -ujQibix) p3<Ji - ^eofci:mWjmfc(a;)crfe , (2-10) 

where the last two terms can be combined into ftb{x) — (l/4)aJcd&(a^)piQ^'^PiQ^'* • It is instructive to see the explicit 
matrix form of the spin connection , 

_f-ieAb \ f -igHb \ 1 fn \ i fn 

— \ n „•„/! + n ^ Ti^Oib n ^ I ^ -r^imk^imb 



-ieAb J \ ig^b J 2 \Q -n ) 4-»""=^""'' \ Q tu 

r 



The upper and the lower rows of this matrix are the spin 
connections for the left (not dotted) spinor and for 
the right (dotted) spinor jy^, respectively. The function 
Aa{x) is the same for both spinors; it is immediately rec- 
ognized as (and will be proved to be) the vector potential 
of the gradient-invariant Maxwell field. 

The signs in the second and third terms alternate be- 
tween the left and right spinors. One can easily under- 
stand the difference in the third term by recalling that 
^abcds'^ is the angle at which the local tetrad is rotated 
in the ab plane when it is transported by infinitesimal 
ds'^. The rotation cjoic is the Lorentz boost; it must have 
opposite signs for the left and right spinors. The spin 
connection ps't^aix) also differentiates between the left 
and right spinors and it will become the field that in- 
teracts with the axial current. This is the most general 
form of a spin connection allowed by Lorentz invariance, 
and all its components are equally important. In order 
for the Dirac field to form compact objects, this general 
form must be restricted. This will be done in the next 
paper Q. 

Next, we have to find the transformation properties of 
the twelve remaining elements of the type ^p~^pa<Jbi'- The 
array 

Ja = Ip^psaalp = ^p^{pz,ai)lp = IpJ^^alp 

is known as the axial current. Using Ta of Eq. (|2.1U|) we 
find that 5j7j = ujbcaJ'^ds'^ . Hence, the axial current jTf, 
is indeed transported as a Lorentz (pseudo-)vector. The 
two other densities, S = ip'ip = 4'^Pi4' and V = ijj^^ijj = 
'0"'"P2'0j are expected to behave as Lorentz scalars. This is 
indeed the case in the sense that their parallel transport 
does not depend on the parameters of the Ricci rotations, 

5S = -2gVKds'', and 6V = 2gSKds''. 



However, in the course of parallel transport the axial field 
H mixes them, i.e., 

(Hereafter, this work substantially deviates from Fock's 
paper 0. Fock requires 5S — and SV = though 
there is no observables that support these conditions.) 
The skew-symmetric Lorentz tensor Mab is supposed to 
be transported as 

SMab = UJacdM\ds'^ + LUbcdM^'^ds'^, 

and an explicit calculation shows that the increment of 
this tensor with the components 



and 



has the correct Lorentz part. However, this increment 
also acquires additional terms due to the field psHa, 

6MQi = {uJOjbMjb - UJi,nbMam)ds'' + geijmMjrn^bds^ , 



6M,k ^{-UJ^cbM\ + uJucbM^'')ds^ -2ge,kmMo^^'^bds'' . 

Therefore, contrary to naive expectations, the parallel 
transport of scalar, pseudoscalar and skew-symmetric 
tensor results in additional terms, which are proportional 
to = c^psHq. This field does not affect parallel trans- 
ports of the vector current and axial current (antisym- 
metric tensor of rank three). However, the axial field H 
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does affect parallel transport of the skew-symmetric ten- 
sor Mab of rank two and its scalar invariants S and V , 
all of which are built with the aid of the off-diagonal ma- 
trices pi and p2- These matrices mix the left and right 
spinors in the equation of motion and various obscrv- 
ables. Mathematically, the two additional fields, Aa{x) 
and correspond to extra degrees of freedom that 

are left to spinors by Lorentz invariance of the tensor 
observables. 



C. Parallel transport and the Lorentz 
transformations. 

The redundant invariance that has appeared in the law 
of parallel transport must show up in Lorentz transfor- 
mation of spinors also. Indeed, let us derive the law of 
the Lorentz transformation for a spinor, ■0 — > Sij:, "0^ — > 
Tp'^ , starting from the same current ja = tp'^aatp, 
which is transformed as ja — > Aj'{x)ji,. This gives 

S+aaS = Aa'{x)ab , (2.12) 

which is the equation that determines S. The matrix 
S{x) can have one of two forms, 

" ( (A+V' ) ' = ( A* ) 5 

A= a5-/37=l, (2.13) 

The "dotted" two-component spinor — rj must be 
transformed either as dual to 0^ = ^ (i.e., by means of 
matrix (A"*")"^ = T2X*T2) or as the conjugated to it (by 
means of matrix A*). Only then wiU (1/2,0) and (0, 1/2) 
be the two different representations of the Lorentz group. 
This transformation should preserve Eq. H2.7|) in its co- 
ordinate form. This form can be derived by introducing 
ctfiix) — e'^^aa and = e^Fa, and accounting for the 
fact that the tetrad vectors are covariantly constant, 

D^el = V^e^ - e>^,,e^ = 0, V^e^ = dxe^ - V^^^el 

(which is just another way to define the Ricci coeffi- 
cients). Hence, the matrices a^t(x) are covariantly con- 
stant also , i.e. D^a^ ~ D^ate'l, — abD^e^ = 0, or, 
explicitly, 

Df,a^ = V^a^ + F+a^ + a„F^ = , (2.14) 

which is just a translation of Eq. (|2.7|l from the tetrad 
basis into the coordinate basis. 

The covariance of (|2.14|) implies that the Lorentz trans- 
formation of a spinor, ip{x) — > -0 = S{x)'ip{x), must be 
accompanied by a similar transformation of the the ma- 
trices afi,{x), i.e. aft{x) afj_{x) = S^^{x)afj_{x)S{x). 
Hence, in order that Eq. (|2.14ll be the same in any local 
Lorentz frame, the spin connections should transform as 

Fp -> = Sr^S-^ + (9^5)5-1 , (2.15) 



where S is the solution of Eq. (|2.12() . The first term of 
this formula "rotates" the spinor indices at each point, 
while the second inhomogeneous term adds a gradient. 

To compare S with the law of parallel transport, 
consider the infinitesimal Lorentz transformation, where 
Aab = Sab+'c^ab, and S = 1+T. Then Eq. H2.12|l becomes 

T+tta + aaT = WabUb , (2-16) 

which is exactly Eq. (|2.7() with ojabcds'^ replaced by an 
infinitesimal Wab = —"^ba- The most general solution of 
this equation that has the required form H2.13|l is 

1 i 
T = iea{x) - -wok{x)p3(Tk ~ -eokim'!^rm{x)crk -(2.17) 

One can recognize the same pattern as in Eq. H2.10|l . 
i.e. the left and right spinors are transformed differently. 
Thus, the local Lorentz rotation of the spinor field acts 
in agreement with the transformation of its spin connec- 
tions. Note, however, that Eq. H2.17|l misses the counter- 
part of the term p3'Ra{x)ds°' in the spin connection even 
though the latter is minimally coupled to the spinor field. 
This is a consequence of the above requirement that the 
matrix S has the form of Eq. (|2.13f) . 

III. EQUATIONS OF MOTION AND 
CONSERVATION LAWS 

It is common to begin field-theory calculations from a 
Lagrangian that has an alleged symmetry of the dynam- 
ical system under investigation and to derive the equa- 
tions of motion and conservation laws using a variation 
principle. Because of a new axial field in the covariant 
derivative of the Dirac spinors we have to proceed more 
gradually. Namely, we must establish the equations of 
motion as the first step, then to check their compliance 
to the data and derive the identities that have a form 
of conservation laws, and only after that to write a La- 
grangian that reproduces the same results. 

A. Equation of motion (Dirac equation revisited). 

In view of a new spin connection, the equation of mo- 
tion of the Dirac field has to be revisited. In fact, we have 
nothing else at our disposal that can be used to create 
an equation of motion, except for the covariant derivative 
of the spinor field. In a tetrad basis it was found to be 
Dip = Daipds"- ~ {daip — Fa'(/')(is° , with the connection 
given by Eq. (|2.1(J|I . The structure of this spin connection 
is two-fold; its spin indices are being used to parameterize 
rotations of the local tetrad basis by means of Pauli ma- 
trices. Its Lorentz index shows the direction of parallel 
transport. The first step is to parameterize the Lorentz 
index a by a spinor and thus to convert this derivative en- 
tirely into a spinor representation. This should be done 
exactly as for the covariant coordinate vector Xa- Ac- 
tually, we have to use Xa x^g = Ixq + f ■ x and 
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Xa — > Xaf3 = Ixq — f ■ X for the left and right spinors, 
respectively. This simultaneous conversion for two com- 
ponents of the Dirac spinor is carried out by means of 
matrix = (1, p^ai) ~ the left and right spinors are 
Lorentz transformed differently. An obvious conjecture 
is that the two spinors do not evolve independently. Oth- 
erwise an electron would be composed of two dynamically 
disconnected fields each having its own light-like current. 
Thus, the equation of motion for the spinor field must 
be compiled according to the correspondence principle, 
which can take different forms. The Dirac form that ap- 
peals to the relativistic = E"^ — = m? can be found 
in almost any textbook. Another form can be traced back 
to the Fock proper-time argument 0- It employs virtu- 
ally the same classical limit but in dimensionless form, 
v? = — vP = \ . (On can even view it as a prototype of 
the commutation relations, a°'pia^ + a^pia"' — 2pig°'^ , 
that sets the Clifford algebra of the Dirac matrices.) It 
is equivalent to the previous one as long as P'^ = mu'', 
but it does not allow for the continuous limit of w? 0. 
This second form is more restrictive - it is a requirement 
that the electron has a rest frame, and a continuous limit 
of ^ just does not exist. Assuming this form, and 
taking the linear relation, w^P'^ = m as a prototype, we 
can write the following version of the Dirac equation, 

a'^(aa^-^a^/')+^mpl?^ = , (3.1) 

with the spin connection again given by Eq. H2.10|l . This 
equation includes an additional vector field Ha(a;). The 
existence of such a field must be confirmed by data at 
the same level of confidence as in the case of the Maxwell 
field. It turns out that the axial field is responsible for a 
dynamic auto-localization of the Dirac field, so that the 
new equation remains very close to its classical prototype. 
The conjugated equation reads as 

{dai)"^ ~ V'+r+)a" - imtlj+pi = . (3.2) 

The differential operators of these equations are Hermi- 
tian (symmetric) , which is confirmed by the conservation 
of the vector current . The question if they can be self- 
adjoint in the presence of an axial field is not trivial and 
will be addressed later on. 



B. Conservation identities for vector currents. 

The equations of motion 1)3. l|l and H3.2|) immediately 
lead to two well-known identities. 



- ^ 0, (3.3) 

and 

V^J^ = -La^lV^^+Paa^'V-] = '^mV. (3.4) 



The first of them, indicates the conservation of the prob- 
ability current. The second one states that the axial 
current cannot be conserved, in a full compliance with 
its space-like nature. A major factor here is a finite 
pseudoscalar density (the mass m is no more than a 
scale parameter), which is a measure of a dynamic inter- 
play between the left and right components of the Dirac 
spinor. This factor become active when a Dirac particle 
experiences an acceleration and its field deviates from a 
perfect parity-even configuration encoded in the conven- 
tional Dirac equation (Cf. |1,0)- 

C. Energy- momentum conservation. 

Physically, the energy momentum tensor is supposed 
to be designed as the flux of the momentum density (a"" 
is the quantum mechanical operator of the velocity, and 
the operator Z?^ has the kinetic momentum P^ = nnu^ 
as the prototype). For now, we begin with a commonly 
used form 

n = \ (^^""^^ - ^'■^«"^) • (3-5) 

which looks like a real quantity and is often believed to 
be self-adjoint. To avoid a premature discussion of the 
complicated issue of self-adjointness that will arise later, 
I shall not dispute this opinion, though some of the re- 
sults will indicate that, in general, this belief is wrong 
because the axial potential in D^j^ can be singular. Then 
it is not safe to have a quantum mechanical operator that 
simultaneously acts in two adjoint spaces without taking 
care of a self-adjoint extension of this (formally symmet- 
ric) operator. This will be done in the next paper Q, 
where the condition of self-adjointness will be derived. 

By virtue of the equations of motion, the following 
identity holds 

+ ^V,V^(V'+a"V') • (3.6) 

The last term in this equation is the result of the follow- 
ing transformation of an anti-Hermitian form, 

^P+a^D^i; + i^+Dfa^il) = D a" ip) = V^(V'+a''?A), 
where a purely geometric relation (|2.14|) . D^a^ = 0, is 
employed and neither the equations of motion nor a spe- 
cific form of the spin connections were used. The com- 
mutator in the brackets is the curvature tensor, 

= a^r, -9,r,, + r,r^-r^r,. (3.7) 

Since the spin connection Fa is defined in terms of its 
tetrad components, the primary objects are the tetrad 
components of this tensor. 

Dab = eaDcpe", = DaD^ - DbDa + C'ab^c , (3.8) 
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where, 

DaDb - DbDa = dbTa - daTb + TaTb - TbTa 

and the structure constants are C^^ = w'^^j, — to'^ba- '^^'^ 
commutator of covariant derivatives can be expressed 
in terms of two gradient invariant tensors (the field 
strengths) , 

Ff,, = d^A, - d,A^, U^, - a^N, - (3.9) 

and the Riemann curvature tensor. Since the vector cur- 
rent is conserved, we have 

where Ra^i = R^a X^i Ricci tensor. The final answer 

reads as follows. 

The last two terms on the right hand side are imaginary 
and (if written in natural units) proportional to Planck 
constant. These terms exactly cancel each other. 

So far, we have the divergence of the spinor energy- 
momentum tensor on the left and the gradient-invariant 
Lorentz forces from the electromagnetic field and a new 
axial field on the right. The Lorentz force eFg-^j'' is 
due to the conserved probability current ja and it drives 
the left and right currents in one common direction. The 
axial force, gUcr^J'^, is due to the not conserved axial 
current J7q, and it pulls left and right currents in oppo- 
site directions. This new force must act between any two 
Dirac particles. Its effect is parity-odd; it has not been in- 
troduced ad hoc and it cannot be rejected a priori within 
a realm of various PNC atomic or nuclear phenomena. 
The extra term which is proportional to the axial poten- 
tial Hp is not gradient invariant. 

The next step is to convert this equation into the di- 
vergence of one common energy momentum tensor for 
all fields in the system. Hence, we need the equations 
of motions for the fields Afj_ and H^. Eqs. 1)3. 9|l immedi- 
ately yield the first (without the sources) couple of the 
Maxwell equations for the field strengths F^,y and U^,y, 

VAe'^'^^Fp. = , VAe"^^"Z^p. = . (3.11) 

The equations that interconnect the fields and currents 
must be conjectured from the existing data. We may 
rely on the following facts: The hydrogen spectrum in- 
dicates that with great accuracy the electron moves in 
the field with potential U — —e/r that satisfies the equa- 
tion AU = — 47re(5(r). We can also refer to the data on 
PNC phenomena as an indication that the axial current 
is involved in electron dynamics. 



The potentials and are Lorentz vectors. There- 
fore, if we identify U — Aq and — ei5(r)=jo in the rest 
frame of the atom, then the only possible Lorentz invari- 
ant form of the second couple of the Maxwell equations 
is 

V^F'^^ = ej'" = eiV'+a'V] , (3.12) 

so that F^^ is the massless gradient-invariant Maxwell 
field. Actually, this choice is straightforward and simple 
because the current is conserved. Then, we can present 
the gradient-invariant Lorentz 4-force of the electromag- 
netic field as the divergence of the energy-momentum ten- 
sor, 

qA ^ ^ px.p^^ ^ {l/4)6^FP''Fp, . (3.13) 

For the axial field H^, the choice of the second equation 
is far less obvious because, except for the "Lorentz force" 
gl^a^S'^ , there is an additional not gradient-invariant 
term ^grrii^^ij}^ p2'il} = ff^Ai(Vi/ J'"'). The data on atomic 
spectra, which indicate that the external field does 
not depend on time and does not deviate from the spher- 
ical symmetry, facilitates this choice. The only non- 
vanishing components can be Ho(?') and H,.(r). Then, 
hii^r = —dr'i^o{r) is the only non-zero component of the 
field strength. This field must be strongly confined near 
the nucleus, which requires a scale. A plausible form is a 
Yukawa potential, so that 

-V^Mr) + M^Mr) ^ gJoir) , 

which, by the argument of Lorentz invariance, suggests 
that Hp is a massive neutral vector field, 

VM"^ + M^H^ = gji" = gitp+paa^ij] . (3.14) 

(We have to put the mass term here to avoid a confiict 
with the anti-symmetry of U^^. Indeed, S^J''' ^ 0.) 
Taking the covariant derivative of this equation we get 

M^VpH^ = g\7^,J^' ^ 2 g mV . (3.15) 

Using (|3.14|) and (|3.4I) one can transform the Lorentz 
force of the axial field in (|3.10l) into the divergence of its 
energy-momentum tensor, 

5[?/.+P3a"V]W.M + 2g"^Hp V+P2^ = ^xO\ , (3.16) 
where 

+S^^(^^UP-^U,, + \M^^''^py (3.17) 
Putting Eqs. H3.10|l - H3.17|l together, one finds that 

VA(r\ + e\ + e\) = o, (3.18) 
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i.e., the total energy-momentum of three interacting 
fields, -0) ^/Ji and is conserved, which is an addi- 
tional indication that the system of equations of motion 
is self-consistent. In the next section we shall address 
the problem of the hydrogen atom and find the effect 
of the static fields Ho(r) and H, (r) on the solutions of 
the Dirac equation. Of particular interest will be the 
quantum numbers of this problem and the existence of 
stationary states. 



IV. HYDROGEN WITH THE DIRAC PROTON 

An obvious way to test the existence and properties 
of the new axial field is to study its effect on atomic 
electrons treating the field H^(a;) as a classical external 
field. We shall look at the nucleus as a cluster of Dirac 
spinors. The electric and axial currents of the nucleons 
are the sources of electric and axial fields. The axial field, 
which has been introduced on quite different premises, 
very much resembles the neutral gauge field of the stan- 
dard model of electro-weak interactions. The neutron- 
rich nuclei are known to have a larger weak charge. For 
the sake of simplicity, the field of the nucleus will be taken 
as time independent and spherically symmetric. The hy- 
drogen atom is used as an example where a new pattern 
naturally emerge in the course of solving the Dirac equa- 
tion with axial the field H. 

For hydrogen-like atoms we expect a perfectly spheri- 
cally symmetric geometry with a natural metric 

ds^ ^ g^^dx^'dx" = dt^ - r^d0^ ~ sin^ 0d(p^ - dr^ . 

The tetrad vectors form a diagonal matrix, and the only 
non-vanishing components of the Ricci rotation coeffi- 
cients are 



1 COS0 
Wi22 = — ; — 7- 



_ 1 
r 



tJ311 = --. (4.1) 
r 



Assume that in addition to the Coulomb field Ao{r), the 
nucleus is also a source of a spherically symmetric axial 
field with the temporal component Ho(r) and a radial 
one ^riT). These are the only components that can be 
present in a spherically symmetric static object. The 
Dirac equation is 

\ido - eAo - api^a + goi^r ~ ip-i<^z{pr + -) 



\coU 



i-^^d^ - mpAi^ = Q . (4.2) 
r sm 



In terms of a new unknown function, ip{r,9,ip) 
rVsinOip , this equation becomes 

[ido - cAq - g(p3^o - 0-3^^) -I- ps.(7z{-idr) 
0-2 



H {-laide 

r 



sin6' 



d^)-mpM = Q. (4.3) 



treats spinor field according to its original definition as a 
Lorentz spinor and a coordinate scalar. (A conventional 
form, with spinors related to Cartesian coordinates, fol- 
lows from this one, if at each space-time point the spinor 
is Lorentz- rotated at a corresponding angle.) In order 
to find its solution one has to separate the angular and 
radial variables in this equation. This is known to be 
a somewhat tricky problem even in the pure QED case 
(when Ha = ). 

The Hermitian operators in Eq. (|4.3|) are the tetrad 
components of the momenta = —idr, pi — —ir~^dg 
P2 = —i{r sin 9)~^ dip. The operators pi andp2 are clearly 
associated with the angular motion. If the coefficients 
in this equation where not the matrices, it would have 
already been an equation with the variables separated, 
which would match the perfect spherical symmetry of 
the external fields and H^. The problem is that the 
operators of radial and angular momenta do not commute 
(they anti-commute, [a^ps, (aipi -I- a2P2)]-i- = ). A 
regular way to avoid this obstacle is as follows 0, |^ • One 
attempts to construct a minimal set of the operators that 
commute with the Hamiltonian. For example, one can 
check that the commutator [013^3, pi{aipi + a2P2)]- = 0, 
and take the operator pi{aipi + 02^2) as a generator of 
the conserved quantum number. This trick works when 
H = and it is very instructive to see the details of its 
failure when H ^ 0. 

The conventional operator of the angular momentum 
is £ = ['rxp\+a/2. An additional operator C — (?•£— 1/2 
commutes with the orbital momentum, [C, {r x p)] ~ 0, 
and has the properties, C{C — 1) = [f x p]^ and ~ 
+ 1/4. Therefore, if k is an eigenvalue of operator £ 
we obviously have k{k — 1) = l{l + 1) and > 0. On 
the other hand, if C = piC, then (£')^ = and these 
operators have the same sets of eigenvalues. In the tetrad 
basis, these operators are 



^^5^): C,^-^^^ + ^a^ . (4.4) 



Equation H4.2|) is the Dirac equation in the tetrad basis. 
This is a primary form of the spinor equation because it 



In terms of the auxiliary operator £' (which has the same 
set of quantum numbers as the operator of the angular 
momentumhut is a different operator) and the projection 
£3 of angular momentum, the Dirac equation becomes 

[ido - eAo - 5(^3^0 - cr3H,.) -I- P3.(7s.{-~idr) 

C 

-P20-3 mpi]%jj = 0. (4.5) 

r 

If = then the operator C commutes with the Hamil- 
tonian and we can require the wave function to be an 
eigenfunction of the Hamiltonian and these two opera- 
tors, 

C'tp = Ki), and Csip ^ {m^ + 1/2)^. (4.6) 

This is a complete solution to the problem. The system 
has two "angular" quantum numbers and a radial one. 

Because the presence of the components Ho(r) and 
Hr(r) preserves the spherical symmetry, we can try to 
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look for the general solution of the following form (the 
notation will become clear later on), 



,(4.7) 



The properties of this form are easy to analyze after we 
substitute it into Eqs. (|4.5|l and (|4.t)|) and put the results 
together (in each line the first equation {—) belongs to 
(|4.5|) and the second equation (==) is taken from (|4.6|) '). 



{[tdo - eAo - <?(Ho + N,) - tdrUL^y - mhn^y + 
{[ido - eAo - .g(Ho - H^) + idr\hLiZ = mfniZ 4 
{[ido - eAo + 5(^0 - H,) + idr]hR^y = mfL^y - 
{[ido - eAo + g{^o + K) - i9r]/i?j.-Z = mhLiZ - 



-{do J-Td^)hLiZ == {ni- i-)hR^y , (a) 

r smO r 

-ide + -^d^)fL^y 
r smu 



(m + *-)/flx2, (b) 



[de - —,d^)fRiZ == (to + i-)h\y , (c) 
smfc' T 

{de + -l^d^)hu^y == (m - i-)hLiZ , {d) 



smf 



(4.8) 



One can easily see that these equations do not go well 
together. Indeed, after the spinor of the form (|4.7|l is 
substituted into Eas. (|4.6|l it becomes apparent that the 
angular variables can be separated only when /l| ~ Jri 
and /i/jf = ft-L^, which is not consistent with the pres- 
ence of the field in Eas. H4.5|l . In this case, the Dirac 
equation with the axial field breaks up into two different 
(and thus, incompatible) systems of equations for only 
two radial functions. 

Nevertheless, just by inspection, one can verify that 
the angular functions yk.m{0, v) ^-nd Zk^m{0, f) that sat- 
isfy the equations 



{de 



-,d^)2k^m{d,'P) = ~kyk,m{0,v), 



{dg + —d^)ykMG, <f) = kZk,UO, f), (4.9) 



do separate angular variables in the Dirac equation (|4.5|l . 
and do not separate them in Eq. 1)4. The angular de- 
pendencies 3^fc.m(6', V?) and Zfc.,„(6', 93) are split between 
the components of Dirac spinor that correspond to the 
opposite (outward and inward) directions of the radial 
component of the electron spin, an eigen-state of the op- 
erator (73 associated with the tetrad vector e^. 



After this separation of angular variables, Eqs 
yield the following equations for the radial functions 



{[ido - eAo - g(Ho + K) - idr\fL] = mhn^ - i-hn , 

k 

{[ido - eAo - 5(^0 - %) + i^rJ/iLi = to/hx + i-^L\ , 

fc 

{{ido - eAo -f 5(^0 - H^) + idr\hR,^ = mfL] + i-fni , 



{[ido - eAo + 5(^0 + "^r) - idr]fRi = ^^^li - i-hR^ 



{a) 

(c) 
id) 



(4.10) 



r 



Now, it becomes clear that the axial field in spin connec- 
tion causes the qualitative change in the theory. Indeed, 
in the absence of the axial field H, we have two identical 
couples of equations for the radial functions that yield the 
well-known energy spectrum of the hydrogen atom. In- 
terestingly enough, this is exactly the spectrum given by 
the Sommerfeld formula of the fine structure, i.e., the so- 
lution of the relativistic Keppler problem with the Bohr- 
Einstein quantization of of the radial motion and of the 
precession of perihelion. ( The Bohr-Sommerfeld addi- 



tional quantum number k corresponds to the quantiza- 
tion of this precession, and within this classical model 
there is no theoretical arguments that would require us 
to reject the linear oscillatory orbit with k = 0. It is 
excluded only because the spectral data provide no ev- 
idence of stationary states with k = in hydrogen-like 
atoms.) Therefore, without the axial field, the quan- 
tum eigenvalue problem yields a classical relativistic an- 



swer, in accordance with the prototype u^P^^ 



of 



the Dirac equation. As long as the spinor polarization is 
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balanced between the left and right components, there is 
not any visible spin effects and the electron moves in an 
intermediate domain were the Maxwell field absolutely 
dominates. The effects of the axial field are genuinely 
small and proportional to the probability to find the elec- 
tron inside the nucleus. The mixing of levels that signals 
the possible parity non-conservation is reasonably well 
described by perturbation theory in the non-relativistic 
limit 

In the presence of the axial field, the spinor polar- 
ization becomes more agile and two new physical pat- 
terns develop. First, the angular functions (|4.9I) are not 
connected in any way with the angular momentum of 
the Dirac field [llj. Therefore, the angles 6 and ip now 
belong to the internal space of parallel transport of the 
spinor field along a closed spherical surface. This, invisi- 
ble from the outside, internal space is very likely to have a 
non-Abelian symmetry group and, possibly, a local gauge 
structure associated with it. Therefore, the well-known 
non-Abelian theories may have a realization within this 
scheme. 

Second, since the angular functions y{9, ip) and Z{9, ip) 
are not connected with the angular momentum there is 
no reason to exclude the case fc = 0, which indicates 
that a new physical pattern develops at short distances. 
This choice corresponds to the physical situation when 
the Dirac field is radially polarized in the internal geom- 
etry of a spherical shell, so that the parallel transport 
of a Dirac spinor along a sphere does not change it. A 
well-known prototype of such a state is a fully occupied 
electron shell of a noble gas, which is symmetric to the 
extent, that none of the electrons can be assigned in- 
dividual quantum numbers associated with the angular 
momentum. The only topologically distinctive direction 
is the radial one. For such states, the Dirac-type equation 
(|4.8() again splits into two different independent systems 
for two polarization modes of the Dirac wave function. 
One of these modes has only two non-zero components, 
^1 and ?7i, which correspond to the positive (outward) 
sign of the radial polarization. The second mode has 
only ^2 and 772 components. Its radial polarization is di- 
rected towards the center (inward). Let us look for the 
stationary solutions of these equations. For the outward 
mode we have Eqs. (|4.1Ul a.c'). 

[E-i- - eAo - gHo - gH^ - idr]fLi = "mhRi, (a) 

[E^ ~ eAo + g^o ~ gK + idr]hR^ = mfL^, (c) (4.11) 

For the inward mode (Eqs. ^A.W\ b.d) ) we have the same 
system with the opposite signs of Hp and H^. 

The wave functions of the inward and outward polar- 
ization modes have different time dependences. It is in- 
structive to see by how much these new solutions deviate 
from a purely Coulomb problem. The equation for the 
static component of the axial potential produced by a 
compact spinor configuration will be derived in the next 
paper 4], and it will yield the expression (|4.12|) below. 
At the moment, we shall take a short cut by noticing 
that in a perfectly spherical static case the field H can 



have only two components Ho(r) and ^^(7"), so that only 
Uor{r) ^ 0. Then by virtue of H3.14|l and I3.15|l we have 
a short-range Yukawa field Ho(r) with the source gJo{r) 
which is spread at a distance ~ ^/M, and a radial 
field Hr(r) for r > rmax, 



gMr) 



Q{r) 



1 2g'^m 
^ M2 



Vir)r'^ dr. (4.12) 



The latter is determined from "Gauss' law" H3.15|l for 
the field H with the nuclear pseudoscalar density as a 
source. Depending on how singular the distribution of 
this density is, the shape of the potential Hj. can range 
from nearly homogeneous within the nucleus to as singu- 
lar as r^2. 

The position of the potential Hq in Eqs. (|4.11() is such 
that one can absorb the Ho(r) into a phase factor by 
means of the substitution, 



{f,h}ti = {f,h}ne^'^^o''"^'^'^ 



(4.13) 



Taking an attractive Coulomb potential with the point- 
like charge + Ze for Ao{r) and Hr(r) from Eg. (|4.12|l 
we obtain two systems of equations for two topologically 
distinct modes, 



ido 



ido 



Ze^ Q(r 
~F ~ 

r 



Ze' 



Q{r) 



+ idr ]h^i = mf^i . 



(4.14) 



Let us assume that Q{r) > Q within some range of radius 
r and take a close look at these equations. The outward 
"f-component of the Dirac spinor is universally attracted 
and the inward J,-component is universally repulsed from 
the positive central pseudoscalar charge. From this per- 
spective, the distribution of the pseudoscalar density in 
the integrand of the "charge" Q is the most intriguing 
issue. If the pseudoscalar density is localized in a very 
small volume, the potential iQ/r^ can overwhelm the 
Coulomb potential. It will be strongly repulsive for the [- 
component and strongly attractive for the |-coniponent, 
even reaching the critical boundary of falling onto the 
center. One may wonder if these dynamics can lead to 
the formation of stable and/or metastable states, and 
what the quantum numbers of these states are and how 
rich the spectrum of these states is. The radial depen- 
dence of gHj.(r) is such that, unless it is screened at very 
short distances, this potential always wins. While the 
Coulomb potential can become catastrophically attrac- 
tive only at Z > 137, when it overwhelms the centrifugal 
barrier, nothing can withstand the attractive —Q/r'^ at 
small r, regardless of the magnitude of Q. ( V. Gribov 
made an attempt to incorporate this phe nomenon into 
his picture of quark/color confinement |1^ . l In this case 
of a singular potential, the conventional wisdom regard- 
ing what is attracted or repelled does not work any more. 
For example, a positive singular potential can be attrac- 
tive for states of negative energy with no lower limited to 
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the energy spectrum. Then the Hamiltonian of the Dirac 
equation is not self-adjoint. This observation will be the 
starting point of the subsequent paper Q , where the ap- 
proach initiated here will be reconsidered once again with 
special attention to the self-adjointness of the Dirac op- 
erators and to the existence of meaningful dynamics in 
the presence of so singular potentials. 



V. NON-RELATIVISTIC LIMIT 

At the atomic scale, the relativistic effects are small 
corrections unless polarization (spin) effects are involved. 
A formal v/c expansion of the Dirac equation naturally 
leads to the modified form of the Pauli equation. Let us 
proceed from the Dirac equation with the axial field (|3.1|) 
in Cartesian coordinates, where the Ricci coefficients are 
absent, 



C Q \ 

ido - —Ao - -7-^3^0 -0 
he he J 

+a(^iV-^I-^P3H^)V' = ^PiV'- (5.1) 



In terms of the large and small spinor components, 

(/) = e~™'='*/'"X^ + 77), and x = e"""'^'*/'"^^ - v) , 
(i.e., in standard representation) these equations become 

[ihdt-eAo- g{a-t)](j) 

= [a-{cp + eA)+g^a]x , (5.2) 
[ihdt - eAo - g{a ■ T?) -I- 2tooc^] X 

= [3 ■(cp + eA)+g-^Q\(i3, (5.3) 

where p = —ihV. In the lowest approximation with re- 
spect to 1/c, equation 15.3|l vields 

X = T7^['?-(p+-^)+5^o]0, 
Zmc c 



and inserting this approximation into Eq. (|5.2I) we arrive 
at the Pauli type equation. 



Ot Zm c Zmc 



igh 
2mc 

eg 



[((7-V)Ho + No(<T- V)]0 
^o(a-A>+^Hg</,. (5.4) 



The first three terms on the right-hand side are the well 
known terms of the Pauli equation. Despite the presence 
of the axial field, the kinetic Schrodinger type term cap- 
tures only the convection pattern and it is not altered by 
the axial field. The fourth and the fifth terms are due to 
the vector and scalar parts of the axial potential, respec- 
tively. The latter is frequently used in non-relativistic 
calculations. It is attributed to the interaction of ax- 
ial electron current with the vector part of the electro- 
weak current and gives a major contribution to the PNC 
atomic phenomena in heavy atoms. The term g{(7 ■ N ) 
corresponds to the interaction between the electron and 
nuclear spin. A detail analysis of spin polarization ef- 
fects that lead to PNC in atomic physics can be found 
in the book If ^ is considered in the wider context 
of an arbitrary external field then one can introduce an 
effective, with respect to its action on atomic systems, 
magnetic equivalent of this field. 



h e fiB 



(5.5) 



Its effect is qualitatively indistinguishable from a true 
magnetic field. 
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